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mpnnc A f Ut0matlC differentiation (AD) is a powerful computational method that provides a 
C0I ^ tin £ f xact sensitivity derivatives (SD) from existing computer P programs 
for use in multidisciplinary design optimization (MDO) or in sensitivity analysis 8 The 
Mathematics and Computer Sciences Division of Argonne National Laboratory and the 
Center for Research on Parallel Computation at Rice University have developed a ore 
compiler AD too for FORTRAN programs called ADIFOR. The ADIFOR too^ has been 
asily and quickly applied by NASA Langley researchers to assess the feasibility and 
computational impact of AD in MDO with several different FORTRAN programs These 
include a state-of-the-art three-dimensional multigrid Navier-Stokes fl 0W P solver for wings 
or aircraft confirmations in transonic turbulent flow. With ADIFOR the user specifies 
sets of independent and dependent variables within an existing computer code ADIFOR 
then .races the dependency path throughout the code, apSZn lto formulate 
derivative expressions, and generates new code to compute the required SD matrix Thp 
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What are Sensitivity 
Derivatives? 

♦ Sensitivity Derivatives (SD) describe how one 
thing changes with respect to another thing 

♦ Example: 

How a car’s speed changes when braking 

- sfowly at first, then more quickly 
(how much ) 

- speed decreases as braking increases 
(which way ) 

♦ £D’s describe how much and which way to 
change the variables in a multidisciplinary 
design optimization (MDO ) 
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Objectives 


♦ Obtain exact SD using the computational 
technique of Automatic Differentiation ( AD ) 


♦ Assess the feasibility and computational 
impact of AD in a typical MDO problem 


The SD Matrix 
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♦ Sample inputs: Mach number or geometry 

♦ Sample outputs: wing pressure coefficients or grid 

♦ £D matrix required in MDO 
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Calculation of the SD Matrix 


♦ Divided differences (DD) (baseline 4- perturbations) 

— Proper step size difficult to determine 
-Truncation & resolution errors 

♦ Hand coding (quasi-analytfcal) / symbolic 
manipulators 

-Manual dependency checking 
- Error prone and time consuming 

♦ Automatic Differentiation (AD) 

— Automatic dependency checking and derivative 

coding 

— Exact derivatives via chain rule 
-Quick & easy; possible speed-up over DD 


Th© AD Tool 


♦ ADIFOR (AD of FORTRAN) / PARASCOPE 
(Argonne National Laboratory and Rice 
University) 

♦ User identifies dependent and independent 
variables in program 

♦ ADIFOR follows program flow, traces 
program dependency paths 

♦ ADIFOR formulates exact derivatives via the 
chain rule 

♦ ADIFOR generates new code for derivative 
objects 
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The CFD Codes 


♦ WTCO: wing C-0 grid generation 

-Algebraic 

-Transfinite interpolation 

♦ TLNS3D: 3-D thin-layer Navier-Stokes solver 

- Finite-volume, central-differencing 
-Grid sequencing, multigrid 
-Scalar artificial dissipation 

- Baldwin-Lomax turbulence model 


ADIFOR Applications in CFD 


♦ WTCO wing grid generation program 

-Independents: thickness, cmax, twist 
- Dependents: grid coordinates (x, y, z) 

♦ TLNS3D Navier-Stokes flow solver 

-Independents: grid coordinates (x, y, z) 
-Dependents: pressure coefficients (C p ) 
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ADIFOR Applications in CFD 


♦ WTCO wing grid generation program 

-Independents: thickness, cmax, twist 
~ Dependents: grid coordinates (x, y, z) 

♦ TLNS3D Navier-’Stokes flow solver 

-Independents: grid coordinates (x, y, z) 

- Dependents: pressure coefficients (C p ) 

♦ WTCO-TLNS3D coupling via file transfer 

-Grid 

-Grid SD matrix 

- Application of chain rule 

d(Sect) d(Grid) d(Sect) 


Computational Results 

♦ ONERA M6 wing planform 

♦ NACA 2412 airfoil sections 

♦ 97 x 25 x 17 grid 

♦ = 0.84, a = 0.00, Re = 1 1.7 x 10 6 

♦ Wing C p and SD of wing C p 

♦ Coloring: white/red = large, blue/black = small 

♦ Several geometries: 

- baseline 

-thickness perturbations ± 

- cmax perturbations ± 

-twist perturbations ± 
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Summary 

♦ Feasibility of using AD in CFD demonstrated 


♦ A PIFOR calculated exact geometric SD for 
grid-flow coupling similar to MDQ problem 


♦ APIFOR calculated SD through complex 
algorithm for nonlinear problem 


♦ APIFOR processing easier & faster than 
quasi-analytic method 


♦ AjD competitive with & more accurate than 
divided differences 
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Sensitivity Derivatives = 
BETTER 

Multidisciplinary Design 
Optimizations: 
PRODUCTS & PROCESSES 

Automatic Differentiation = 
EASILY, QUICKLY & RELIABLY 
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